SQUARE FUNCTION AND HEAT FLOW ESTIMATES ON 

DOMAINS 
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Abstract. The first purpose of this note is to provide a proof of the usual 
square function estimate on LP{Q). It turns out to follow directly from a 
generic Mikhlin multiplier theorem obtained by Alexopoulos, which mostly 
relies on Gaussian bounds on the heat kernel. We also provide a simple proof 
of a weaker version of the square function estimate, which is enough in most 
instances involving dispersive PDEs. Moreover, we obtain, by a relatively 
simple integration by parts, several useful L''{Q; H) bounds for the derivatives 
of the heat flow with values in a given Hilbert space H. 



1. Introduction 

Let be a domain in R", n > 2, with smooth boundary dfl. Let Ad denote the 
Laplace operator on H, with Dirichlet boundary conditions, acting on L^(r2), with 
domain H^{n)nH^{n). 

The first result reads as follows: 

Theorem 1.1. Let f e C°°{n) and * e C^(R*) such that 
(1.1) ^*(2-2^A) = l, AeM. 

Then for all p € (1, oo) we have 

1/2 



(1.2) \\fhHn)^Cj,\\[Y.\^{-2-'=An).f\'y 
where the operator ^{~2^^^ Ajj) is defined by (I3.20p below. 



Lp{n) 



Readers who are familiar with functional spaces'theory will have recognized the 
equivalence Fp'^ « L^, where the Triebel-Lizorkin space is defined using the right 
hand-side of (|1.2p as a norm. In other words, LP{^l) and the Triebel-Lizorkin space 
Fp'^{il,) coincide. Such an equivalence (and much more !) is proven in [25l [26l [27] . 
though one has to reconstruct it from several different sections (functional spaces are 
defined differently, only the inhomogeneous ones are treated, among other things). 
As such, the casual user with mostly a PDE background might find it difficult to 
reconstruct the argument for his own sake without digesting the whole theory. It 
turns out that the proof of (|1.2p follows directly from the classical argument (in R") 
involving Rademacher functions, provided that an appropriate Mikhlin-Hormander 
multiplier theorem is available. We will provide details below. 

A weaker version of Theorem 1 1.1 1 is often used in the context of dispersive PDEs: 
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Theorem 1.2. Let f G C°°{Q), then for dip e [2,oo) we have 
(1-3) <Cp|(^||vI/(-2-^^A,,)/|li,(,,))'^'. 

The second part of the present note aims at giving a self-contained proof of 
()1.3p , with "acceptable" black boxes, namely complex interpolation and spectral 
calculus. In fact, if one accepts to replace the spectral localization by the heat 
flow, the proof can be made entirely self-contained, relying only on integration 
by parts. Our strategy to prove Theorem 11.21 is indeed to reduce matters to an 
estimate involving the heat flow, by proving almost orthogonality between spectral 
projectors and heat flow localization; this only requires basic parabolic estimates 
in LP(il), together with a little help from spectral calculus. 

Remark 1.3. For compact manifolds without boundaries, one may find a direct proof 
of (|1.3[) (with A £) replaced by the Laplace-Beltrami operator) in [7] , which proceeds 
by reduction to the M" case using standard pseudo-differential calculus. Our ele- 
mentary approach provides an alternative direct proof. However, the true square 
function bound (|1.2p holds on such manifolds, as one has a Mikhlin-Hormander 
theorem from PT| . 

Remark 1.4. One can also adapt all proofs to the case of Neumann boundary 
conditions, provided special care is taken of the zero frequency (note that on an 
exterior domain, a decay condition at infinity solves the issue). The Gaussian bound 
which is required later holds in the Neumann case, see [lOl [9]. 

Remark 1.5. As mentioned before, Theorem 11.21 is useful, among other things, 
when dealing with estimates for wave or dispersive evolution equations. For 
such equations, one naturally considers initial data in Sobolev spaces, and spectral 
localization conveniently reduces matters to data in L^, and helps with finite speed 
of propagation arguments. One however wants to sum eventually over all frequen- 
cies in if possible without loss. Recent examples on domains may be found in 
[IS] or [10], as weU as in [TT]. 

We now state estimates involving directly the heat flow, which will be proved 
by direct arguments. It should be noted that for nonlinear applications, it is quite 
convenient to have bounds on derivatives of spectral multipliers, and such bounds 
do not follow immediately from the multiplier theorem from [2]. We consider the 
linear heat equation on Q with Dirichlet boundary conditions and initial data / 

(1.4) dtu- Adu^O, on fix R+; u\t=o = f & C°° (fl); - 0. 

We denote the solution u{t,x) = S{t)f{x), where we set S{t) — e*^". For the 
sake of simplicity A^i has constant coefficients, but the same method applies in 
the case when the coefficients belong to a bounded set of C°° and the principal 
part is uniformly elliptic (one may lower the regularity requirements on both the 
coefficients and the boundary, and a nice feature of the proofs which follow is that 
counting derivatives is relatively straightforward). 

Let us define two operators which are suitable heat flow versions of 5'(— 2^^-' A^i): 



(1.5) 



Qt = ViVSit) and Q, = tdtS{t). 
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Theorem 1.6. Let 1 < p < +oo, then we have 




Lp(n)j 



which implies, for p G [2, +cxd), 



(1.7) 




and Qt may be replaced by Qt in both statements. 

Notice that there is no difficulty to define Qtf or Qt/ as distributional derivatives 
for / e LP{il), while simply defining '^{—2^^^ Ad) on LP{^1) is already a non trivial 
task. The purpose of the next Proposition is to prove that both operators are in 
fact bounded on LP{fl). 

Proposition 1.7. Let 1 < p < +oo. The operators Qt, Qt are bounded on LP{fl), 
uniformly in i > 0. Moreover Qt is bounded on L^(il) and L°°(il). 

For practical applications, one may need a vector valued version of Proposition 
11.81 Let us consider now u — (M;);G{i,..,Af} for N > 2, where each ui solves ()1.4|) 
with Dirichlet condition and initial data // . Let H be the Hilbert space with norm 
WuWjj = J2i ^iid L'P{Q.;H) the Hilbert valued Lebesgue space. Then we have 

Proposition 1.8. Let 1 < p < +oo. The operators Qt, Qt are bounded on 
LP{yi;H), uniformly in < > and N. Moreover Qt is bounded on L^{il;H) and 



Remark 1.9. One may therefore extend the finite dimensional case to any separable 
Hilbert space. The typical setting would be to consider the solution u to the heat 
equation with initial data f{x, 9) ^ Lg = H . Notice that the Hilbert valued bound 
does not follow from the previous scalar bound; however the argument is essentially 
the same, replacing | • | norms by Hilbert norms. 

Remark 1.10. A straightforward consequence of Propositions 11.81 and 11.71 is that 
the Riesz transforms Ad)"^ are continuous on Besov spaces defined by the 
RHS of (|1.7p : these spaces are equivalent to the ones defined by the RHS of l|1.3p . 
see Remark l3 . 5 1 later on. 

Alternatively, one can derive all the (scalar, at least) results on the heat flow 
from adapting to the domain case the theory which ultimately led to the proof 
of the Kato conjecture ([S1[S]). Such a possible development is pointed out by P. 
Auscher in [4] (chap. 7, p. 66) and was originally our starting point; eventually we 
were led to the elementary approach we present here, but we provide a sketch of 
an alternate proof in the next remark, which was kindly outlined to us by Pascal 
Auscher. 

Remark 1.11. The main drawback from p.7|) is the presence of VS{t) on the right 
hand-side: one is leaving the functional calculus of A^i, and in fact for domains 
with Lipschitz boundaries the operator VS'(t) may not even be bounded. As such, 
a suitable alternative is to replace VS{t) by \/dtS{t). Then the square function 
estimate may be obtained following U| as follows: 

• prove that the associated square function in time is bounded by the 
norm, for all 1 < p < 2, essentially following step 3 in chapter 6, page 55 
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in [4]. This requires very little on the semi-group, and Gaussian bounds 
on S{t) and dtS{t) ([llj) are more than enough to apply the weak (1, 1) 
criterion from 4 (Theorem 1.1, chapter 1). Moreover, the argument can 
be extended to domains with Lipschitz boundaries, assuming the Laplacian 
is defined through the associated Dirichlet form; 

• by duality, we get the square function bound for p > 2 (step 5, page 56 in 

i); 

• from now on one proceeds as in the remaining part of our paper to obtain 
the bound with spectral localization, and almost orthogonality (j3.7p is even 
easier because we stay in the functional calculus. One has, however, to be 
careful if one is willing to extend this last step to Lipschitz boundaries, as 
this would most likely require additional estimates on the resolvent to deal 
with the Aj. 

2. From a Mikhlin multiplier theorem to the square function 

The following "Fourier multiplier" theorem is obtained in [2] under very weak 
hypothesis on the underlying manifold (see also [3| for a specific application to 
Markov chains, and [23] for a version closer to the sharp Hormander's multiplier 
theorem, under suitable additional hypothesis, all of which are verified on domains). 
For m S L°°(R"'"), one usually defines the operator m(— A/^) on L^(ri) through the 
spectral measure dE\: 

r+oo 

(2.1) m(-Ai5)= / dEx, 

Jo 

and m{~A]j) is bounded on L^. 

Remark 2.1. One may alternatively use the Dynkin-Helffer-Sjostrand formula as in 
the Appendix, and both definitions are known to coincide on L'^{fl). However, the 
Dynkin-Helffer-Sjostrand formula seems to be restricted to defining to(— A^i) for 
functions m which exhibit slightly more decay than required in the next theorem, 
at least if one proceeds as exposed in the Appendix. 

Theorem 2.2 ([2J). Let m e C^(R+), N e N and N > n/2 + 1, such that 

(2.2) sup \^d^m{^)\ < +O0. 

i,k<N 

Then the operator defined by (|2.ip extends to a continuous operator on LP{Q), and 
sends L^(fl) to weak L^{fl). 

In order to use the argument of [2], we need the Gaussian upper bound on the 
heat kernel, which is provided in our case by [10]. Once we have Theorem 12.21 all 
we need to do to prove Theorem 11.11 is to follow Stein's classical proof from [22^ 
and we recall it briefly for the convenience of the reader. Let us introduce the 
Rademacher functions, which are defined as follows: 

• the function ro(t) is defined by ro{t) — 1 on [0,1/2] and ro(t) = —1 on 
(1/2, 1), and then extended to R by periodicity; 

• for m e N \ {0}, r„(t) = ro(2'"t). 



we thank Hart Smith for bringing this to our attention 
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Their importance is outlined by the following inequalities (see the Appendix in 

[22]), 

(2.3) Cp||^a„r„(i)||LP <(5]|a,„p)5 < ^ a„r„(i)||iP. 

m mm 

Now, define 

where was defined in the introduction. A straightforward computation proves 
that the bound holds for m=t(i,^). Therefore, 

||m±(i,-A^)/|U.(0) < ||/|U.(n); 
integrating in time over [0, 1], exchanging space and time norms, and using (|2.3p . 

±oo 
3=0 

This proves one side of the equivalence in (|1.2p : the other side follows from duality, 
once we see the above estimate as an estimate from ^^(il) to LP{n; P), which maps 
/ to (vI/(-2-2^"Az3)/),ez. 

3. Heat flow estimates 
In order to prove Proposition II . 61 we need the following lemma. 
Lemma 3.1. For all 1 < p < +oo, we have 

(3.1) \\S{t)f\\L.in)^t^ooO, 

(3.2) sup\\S{t)f\\L.in)<\\f\\Lnn)- 
Moreover, 

(3.3) \\snp\S{t)f\U,^a)<\\f\\Lp(n) < WfU^in)- 

Proof: The estimate (|3.2p clearly follows from (|3.3p . which in turn is a direct 
consequence of the Gaussian nature of the Dirichlet heat kernel, see [10]. The same 
Gaussian estimate implies p.ip . However we do not need such a strong fact to 
prove p.2p . which will follow from the next computation as well (see (|3.4p ) when 
1 < p < +00. Estimate (|3.ip can also be obtained through elementary arguments. 
We defer such a proof to the end of the section. 

3.1. Proof of Theorem 11.61 li p ~ 2 the proof is nothing more than the energy 
inequality, combined with (j3.ip . In fact, for p = 2, we have equality in p.6p with 
C2 = 2. We now take p = 2m where to > 2. Multiplying equation (|1.4p by 
and taking the integral over fl and [0,r], T > yields, taking advantage of the 
Dirichlet boundary condition, 



(3.4) I rd,\\u\\l dt+ [ 
P Jo Jo 



T 

\S/u\^\u\P-^dxdt 

n 

T 
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from which we can estimate either \\u\f^pf^^-^{T) < ||/||^p(q) (which is (13. 2p ) or 

< \\u\\l,^^^iT)+p{p-l) f j \Vu\M'"^dxdt. 

Jo Jo 

Letting T go to infinity and using (|3.ip from Lemma |3. II and Holder inequality we 
find 

n pOO P 2. P P — 2 

ii/iiL(o)<p(p-i)(yjy^ \vu\^dtydxy[j^{snp\ur')^^dx)' . 

The proof follows using again Lemma |3. 11 as 

p-2 



IP 



<Cp\\( / |Vu|^dt) ' ||lp(o) ( II sup|m|||lp(o; 
^Jo ' \ t>o 



Note that we may prove the weaker part, (jl.7p . without assuming the maximal in 
time bound, by reversing the order of integration in our argument. This would keep 
the argument for heat square functions essentially self-contained, without any need 
for Gaussian bounds on the heat kernel. 

Remark 3.2. We do not claim novelty here: our argument follows closely (a dual 
version of) the proof of a classical square function bound for the Poisson kernel in 
the whole space, see [22^. 

We have proved one side of the equivalence in (jl.6p involving the Qt square 
function, in the range 2 < p < +oo; we now prove the other side, by duality. Let 
e L'?(ri), with Xjq = 1 — 2/p, and consider 

r r+oo 

I = { \\7u\^dt)(t){x)dx. 
Jn Jo 

Without any loss of generality, we may assume (p > 0. On the other hand, let 
V — |VmP, then 

dtv - Av = -2\V^u\^, 

and one checks easily that dnV = on dfl. Let Sn{t) be the solution to the heat 
equation on Q with Neumann boundary condition, by comparing v and 5„(t/2)w(f/2) 
(formally, take the difference, multiply by the positive part and integrate by parts) 
we have 

0<v< Sn{t/2)v{t/2) = 5„(i/2)|Vii(i/2)|2, 

and therefore 

(3.5) /< 



n n+oo 

2 \Vu\^ Snit)(t)dxdt. 

Jn Jo 



Now, we also have 

dtu^ ~ Adu^ = -2|Vu|^ 

and therefore 

/<-/ / [dt- A){u^)Sn{t)<pdxdt. 
Jn Jo 

From {dt - A){u^ Snit)<j)) = -2V{u^) ■ V5„(i)(/), we get 



/ < 



4sup\u\{ |Q,^|2 )2( / \VSnim''dt) 
t Jo t Jo 
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The bound we already proved with Qt can easily be reproduced with S{t) replaced 
by Sn{t), and therefore, provided q <2, we may use the dual bound on the square 
function of € L'^{^1) and conclude by Holder, using (|3.3p on the first factor. The 
condition on q translates into p > 4, and the remaining 2 < p < 4 are handled by 
interpolation. 

Remark 3.3. Actually, we may directly bound Sn{t) by a Gaussian in (|3.5p . extend 
by outside fl, and use the heat square function bounds in R". This provides a 
direct argument, irrespective of the value of p. 

It remains to prove the equivalence between the Qt square function and the Qt 
square function. For this, we repeat the duality argument but we replace jVup 
by t\dtu\'^. Notice that dtu is also a solution to the heat equation with Dirichlet 
boundary condition, and if u) = |9tMp, 

{dt-A)w^~2\Vdtu\^. 

Therefore, comparing w and S{t/2)w{t/2), 

< \dtu\^ < S{t/2)\dtu{t/2)\^, 



and 



Now, 



/• r + co !• p+oc 

J= / \dtu\h(t) dxdt < 2 / \dtu\hS{t)(j)dxdt. 
Jn Jo Jn Jo 



J<2 / tdtuAuS{t)(j)dxdt 
Jn Jo 

< / tdt\Vu\'^S{t)4>dxdt -2 / tdtuVuVS{t)(t)dxdt 
Jn Jo Jn Jo 

f f^°° f r^°° dt 

< \Vu\^{l+tdt)S{t)(t>dxdt + 2 / \QtuQtuQt(l)\—dx 
JnJo JnJo i 

from which we can easily conclude by Holder (using Lemma fLSl to bound tdtS{t)4)). 
Duality takes care of the reverse bound, and this concludes the proof of Theorem 
II. 6( except for the equivalence between the Qt and Qt Besov norms in ()1.7p : we 
defer this to the end of the next subsection. 

Notice that, at this point, we proved Theorem II. 2[ but with the 5* operator 
replaced by the gradient heat kernel and the discrete parameter 2^^-' by the con- 
tinuous parameter t. The rest of this section is devoted to proving the equivalence 
between the Besov norms which are defined by the heat kernel or the spectral 
localization. 

Lemma 3.4. Let \ < p < +oo. We have the following equivalence between dyadic 
and continuous versions of the Besov norm: 

4EllQ2-/llL(n)< / IIQ*/llL(n)7<3EllQ2-/lli.(n)- 

This follows at once from factoring the semi-group: for 2~^^ < t < 2~^^^~^\ 
write S(t) = S{t - 2'"^^ ) S (2-"^^ ) and use jS^l)- We now turn to the direct proof 
of Theorem 11.21 from the heat flow version. Let 5* € C(j"(R*) satisfying (jl.ip and 
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denote Ajf =^ ^{2''^^ Ad)/, where *(2-2JA^)/ is given by the Dynkin-Helffer- 
Sjostrand formula (see the Appendix, (|3.20p ). From Proposition 11.61 and Lemma 
13.41 we have 

1/2 

(3.6) ll/IU.ai)<3C,(^|l 

and we will show that l|3.6p implies (|1.2p : it suffices to prove the following almost 
orthogonality property between localization operators Aj and Q2-2k: 

(3.7) Vfc,jeZ, ||Q2-..A,/|Up(n) <2-IJ""'=I||A,/|Up(o). 
Then, from {2^^^^''^)k e and (|| Aj/||ip(Q))j e we estimate 

(3-8) X! ll'32-2'»/|lip(a) ^ II Q2-^fc^37llip(n) 

as an * P convolution and conclude using Lemma [TT51 It remains to show p.7p : 

• for k < j we write 

Q2--Aj/ = 2^/22-20-'=) (2-(2fe+i)/2v^(2-(2fc+i))) 

|^2-(2fe+i)A^5'(2-(2fe+i))^^(_2-2JA^)^-(_2-2JA^)/, 

where we set *(A) =^ j^(A), and * G C;f, * = 1 on supp^*. By Lemma 
[TSl the operators Qa-c^-^+D = 2^(2fc+i)/2vS'(2-(2fc+i)) and Q2-(2;=+i) = 
2-(2fe+i)^^5.(2-(2fe+i)) are bounded on LP(f7) and we obtain jSJ]) using 
Corollary [3ll2] for 

• for k > j we set ^i(C) = ^(C)6xp(0, 5'2(C) — ^(C)i ^i^'i use again 
Lemma [3. 151 to write (slightly abusing the notation as 2^^'' — 2^^^ < 0) 

(3.9) S{2-^'' - 2-^^)Ajf = S{2-^'')^i{-2-^^ Ad)-^2{-2-^^ Ad)/. 
Then 

Q2-2>^Ajf = 2-^=-^) (2-^ VS'(2-2J)^ (s'(2-2'= - 2-2J)Aj/), 

and using again Lemma fTTSl we see that the operator 2^^ VS'(2^2^) is bounded 
while the remaining operator (|3.9p is bounded by Corollarv l3.12l This ends 
the proof. 

Remark 3.5. One may prove a similar bound with Q2-2k and Aj reversed, either 
directly or by duality. Hence Besov norms based on Aj or Q2-2k are equivalent. 

3.2. Proof of Proposition [TTSl For Qt, boundedness on all spaces, including 
p = 1, +00, follows once again from a Gaussian upper bound on dtS{t) (see [11] or 
[T3])- However the subsequent Gaussian bound on the gradient \7xS{t) in [TT] is a 
direct consequence of the Li-Yau inequality, which holds only inside convex domains. 
We were unable to find a reference which would provide the desired bound for Qt 
in the context of the exterior domain. Therefore we provide an elementary detailed 
proof for Qt- Furthermore, we only deal with 1 < p < 2 or powers of two, p = 2'", 
rn e N*: complex interpolation takes care of remaining values of p, though one 
could adapt the following argument to generic values p > 2, at the expense of 
lengthier computations. 
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Set v{x,t) — {vi, ..,Vn){x,t) :— Qtf — t^^^\7u{x,t) and assume without loss of 
generality that Vj are real: we multiply the equation satisfied by v by wltil^""^, where 
|wp = J2^=i ■^1' ^'^d integrate over ft, 

(3.10) dJ-\MUn))-J2 [ {{^■y)v,)-v,\vr'da+ 

+ f \Vv\M'-'dx+^^P^ j V{\vnvrUx=Uv\\l 



Lp(n)' 

where ~u is the outgoing unit normal vector to 9f2 and da is the surface measure 
on do,. We claim that the second term in the left hand side vanishes: in fact we 
write 



(3.11) / (^-v^.o-^.br^da 



Ida 

= / 5.(|9.U|2 + \^tangU\^){\duU\^ + tangU\'')'^P-^'^ da, 

^ Jan 

and from wjao = the time and tangential derivative (i9t, Vtang)u|an vanishes; 
furthermore, using the equation, &f,u = on dO,. 

Remark 3.6. Notice that while this term does not vanish with Neumann boundary 
conditions, it will be a lower order term (like |Vup on dVl) which can be controled 
by the trace theorem. 

Now, if 1 < p < 2, multiply by integrate over [0, T], 

df 

\\v\\lnn){T)< \\Qtf\\Un)j<\\f\\l 
where the last inequality is the dual of (|1.6p . Hence we are done with 1 < p < 2. 

Remark 3.7. We ignored the issue of v vanishing in the third term in (j3.10p . This 
is easily fixed by replacing Iwl^*"^ by (vs+K'F)^"^ and proceeding with the exact 
same computation. Then let e go to after dropping the positive term on the left 
handside of dSTTO]) . 

Now let p — 2™ with m > 1: we proceed directly by integrating p.lOp over 
[0,T],to get 

(3.12) l\\v\\l (T)+ f j \Vv\M'-^dxdt+ 
P Jq Jn 



2 , 

On the other hand (recall dO]) ! 



(3.13) -\\u\\l,,,,AT) + {p-l) I I iv7.i^u,i.-w.^.-liifiiP 



^■'-"W^)(r) + (p-i) / / |V«|Vr^dxdt = i||/||^,(^). 
e from (f3T2l) . dSTS]) . 

4Un)iT)< I ^Mi.^n)dt^^ [ \\^u\\l,^,,^ dt < \\\f\ 



P 

li p = 2 the estimates are trivial since from p.l2|) . (I3.13|) . 

1,, ,,o /"^ 1 ,, ,,o . 1 r ..o . 1 
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Now, let p > 4:-, for convenience, denote by J the second integral in the left hand-side 
of (|3.12|) (notice that the third integral is bounded from above by J), hence 

j = r [ \^^u\^\\/urhUxdt = r [ {y^\diju\^){y^\dju\^)^t^ dxdt, 

Jo Jn Jo Jn j J ' J 

and set 

(3.14) Ik= f [ \S/u\^''\u\P-^''t''-^ dxdt where 2 < 2fc < p. 

^0 Jn 

For our purposes, it suffices to estimate the right hand-side of (|3.12p . which rewrites 

(3.15) \[ t^^-'\\^AUn)dt=\u. 

Integrate by parts the inner (space) integral in 1^ , the boundary term vanishes and 
collecting terms, 

(3.16) 

Jn {p-2k + 1) 

By Cauchy-Schwarz the integral in the right hand side of (|3.16p is bounded by 



therefore for fc > | + 1 we have 

< (2fc-l) jijh 

{p-2k + iy ^2fc-f-i- 

We aim at controlling Im by J^~^Ii, for some ij > which depends on m (notice 
that when p = 4, which is m = 2, we are already done, using k ~ 2 \). Set 
A: = H - (2^ - 1) with j < to - 2, 

^ (2"' - (2^+1 - 1)) 1 I 

J2'"-i-(2J-l) S (2J + 1 - 1) -"^2'"-i-(23 + i-l)' 

and iterating to — 2 times, we finally control Ip by J^^^I^, which proves that Qt 
is bounded on L^'(ri). 

We now proceed to obtain boundedness of Qt on LP{VI) from the Qt bound; this 
is worse than using the Gaussian properties of its kernel, as the constants blow up 
whenp — > 1, +00. It is, however, quite simple. By duality is bounded on U'{Q)^ 
and 

Q, = tdtS{t) ^ tS{^)/\S{^) ^ 2 J|5(i)V • a/|v^(^) = 2Q|Q. , 



and we are done with Lemma 11.81 

From the previous decomposition, we also obtain 

IIQt/l|LP(0) ^ \\Qtf\\LP{n), 

which implies that any Besov norm defined with Qt is bounded by the corresponding 
norm for Qt. The reverse bound is true as well, though slightly more involved. We 
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provide the proof for completeness. Consider f,h £ C§°{n) and (/, 5) = In -I'd- 
Then 

if, 9) = - {dtS{t)f, h) dt = ~2 {dtSit)f, S{t)h) dt 
JO Jo 

= 2/ {dtS{t)f,dsS{s)h)dtds = A {VS{s)dtS{t)f,VS{s)h)dtds 



t<s 



< 



VS{t)dsS{s)fdt 



\7Sis)h\\p,ds< / ^s\\dsS{s)f\\p\\VS{s)h\\p> ds 



where we used our bound on \/t.WS{t) at fixed t. Then 

(Lh) < f \\QsfUQsh\\p>- 

J s ^ 

from which we are done by Holder. 



3.3. Proof of Proposition [T771 Let us consider now the vector valued case u = 
jv} for N > 2, where each ui solves (jl.4[) with Dirichlet condition and 
initial data //. For the sake of simplicity we consider only real valued uj, and write 

N n n N 

1=1 j=i j=i 1=1 

Notice that n is the spatial dimension and is fixed through the argument: hence 
all constants may depend implicitely on n, while N is the dimension of H. For 
p = 1, +00, the boundedness of Qf follows from the Gaussian character of the time 
derivative heat kernel, which is diagonal on H. 

We proceed with Qt- Multiplying the equation satisfied by ui by in- 
tegrating over and summing up we immediately get (|3.4p . We now proceed to 
obtain bounds for v{x,t) = {vi{x,t))i, where vi{x,t) = t^^^'Vui{x,t). Multiplying 
the equation satisfied by vi by where |wp — tjVitp, summing up over I and 

taking the integral over fl yields 



(3.17) -|kllL(o)(T^) + EE r / \y(.d,uO\'\Vur'dxdt+ 

+ £ I MvumvurH^/' dxdt = £ \\vu\\i^^,, f/'-' dt = i/. 



where |V(9,ziOI' = EtiidljUi)^ IVup = Ei=i E;=i(5jW/)'- Notice again that 
the boundary term vanishes. Denote the last two integrals in the left hand side 
by Ji, J2. Like before, we perform integrations by parts in Ik defined in (|3.14p to 



12 



O.IVANOVICI AND F.PLANCHON 



obtain 

N 



(3.18) 



N 

^^u\P-^''-^ dx. 



[ \Vu\'^''\u\P'^^ dx ^ -J2 I u/Au;|Vu|2('=-i)|u|P-2fc^^. 
Jn Jn 

N 

-(fc-l)V / WiVMiV(|Vwp)|Vu|2('^-2)|M|P-2'=da; 
1=1 

n . N 

i=l •'^^ 1=1 

For fc > I + 1 we estimate the first term in the right hand side of (|3.18p by 

.AT N 

/ (5^w?)'^'E(Ati/)')'/'|VM|2('=-i)|ur2feda; < 
•'^ 1=1 1=1 

(EE / / iv(a,w,)nvur2dx)i/2( / \vuf>^-p-M^^~^''^^d^f^^, 

j^-^ Jn Jn Jn 
and the second term in the right hand side of (|3.18p by 

(fc _!)(/■ V(a,(|V«|2))2|Vur'ldx)l/2(/' |Vw|4'^-f-2|^|2p-4fe+2^^)l/2^ 

Jn j^j^ Jn 
where we used that 

N n N N 

J2ni\7ui\/{\\/u\') < ^(^u?)^/^(^(9,z.0')'/'|a,(|Vun| < |^.||V^.||V(|VK^|. 

1=1 i=l 1 = 1 1=1 

Since the last term in (j3.18p is negative, while the quantity we want to estimate is 

positive we obtain from the last inequalities 

(3.19) 

/ / |Vup'=|zi|f-2feifc-l d:cdt < {Jl'^ + 4/')/2fe_£_i < (Jl + J2)^/'/2fe-f-l- 

Jo Jo 

From now on we proceed exactly like in the scalar case iterating sufficiently many 
times to obtain the desired result, since we control Ip/2 which is the RHS term of 
(I3l8ll using (I3l9l) . 

3.4. A simple argument for (|3.ip . We now return to the first estimate in Lemma 
13.11 while we only deal with p = 2, there is nothing specific to the case in what 
follows. Let X be a smooth cut-off near the boundary dfl. Then t; = (1 — x)u solves 
the heat equation in the whole space, with source term [x, A]u: 

{1 - X)u = So{t){l - x)uo + I Soit- s)[x.A]u{s)ds, 

Jo 

where Sq is the free heat semi-group. We have, taking advantage of the localization 
near the boundary, 

||[X, A]j.|l^^^^^^ < C(x,x')l|Vu|L.(i.) < +00, 
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by the energy inequality (|3.4|) . The integral equation on (1 — x)u features 5*0 for 
which we have trivial Gaussian estimates, and both the homogeneous and inhomo- 
geneous terms are Ct(L^) and go to zero as time goes to +00. On the other hand, 
by Poincare inequality (or Sobolev), 




which ensures that Hx^^lU goes to zero as well at t — +00. 
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Appendix: functional calculus 

We start by recalling the Dynkin-Helffer-Sjostrand formula ([lUIIl]) and refer 
to the appendix of [18j for a nice presentation of the use of almost-analytic exten- 
sions in the context of functional calculus. In what follows we will also rely on 
Davies'presentation ([12]) from which we will use a couple of useful lemma. 

Definition 3.8. (see [HI Lemma A.l]) Let ^ € C|J°(R), possibly complex valued. 
We assume that there exists ^ e C^(C) such that |9^(2:)| < C|Imz| and = 
Then we have (as a bounded operator in i^(51)) 

(3.20) ^i-h^Ao) = I d^{z){z + h^ADy^dzAdz. 

27r Jc 

The next result ensures the existence of ^ in the previous definition ( see [TBI 
Lemma A. 2] and where it is linked with Hadamard's problem of finding a 
smooth function with prescribed derivatives at a given point): 

Lemma 3.9. // f belongs to Cg°{R) there exists ^ G C^(C) such that ^flM = ^' 
and 

(3.21) |a*(z)| < CAr,vi,|Imz|^, VzeC, VNgN. 

Moreover, if ^' belongs to a bounded subset of (M) ( elements of B are supported 
in a given compact subset of M with uniform bounds ), then the mapping S 9 \1/ — > 
^ e C|J°(C) is continuous and Cn,^ can be chosen uniformly w.r.t ^E* G S. 

Remark 3.10. Estimate p.21|) simply means that d'i'{z) vanishes at any order on 
the real axis. Precisely, if z = x + iy 

In particular if (x) = (1 -I- a;^)^/^ then for any given A > 0, a useful example of an 
almost analytic extension of 4" G C^(]R) is given by 

N 
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where r is a non-negative C°° function such that r(s) = 1 if |s| < 1 and t(s) = 
if |s| > 2. For later purposes, we also set 

N „ 
11*11^ li^^ / 

Our next lemma lets us deal with Lebesgue spaces. 
Lemma 3.11. Let z ^M. and \lmz\ < |Rez|, then Ad satisfies 



(3.22 



for 1 < p < +00, with a constant c = c{p) > and a — a{n,p) > — ^\. 

Remark that, for all h G (0, 1], the operator h^Au satisfies (|3.22p with the same 
constants c and a (this is nothing hut scale invariance) . 

For p = 2 the proof of Lemma lS.lll is trivial by multiplying the resolvent equation 
—Adu + zu = f hy u and we get a = 0; however for p ^ 2 it requires a non trivial 
argument which we postpone to the end of this Appendix. 

Corollary 3.12. For iV > a + 1 the integral (|3.20p is norm convergent and V/i g 
(0,1] 

(3.23) |l*(-/i'Az5)|Up(o)^LP(0) <c||*|lAr+i, 

for some constant c independent of h. 

Remark 3.13. Notice how the Mikhlin multiplier condition (|2.2p on does not 
imply boundedness of ||\E'||Ar+i: we need extra decay at infinity. 

Proof: By scale invariance it is enough to prove p.23p for h = 1. The integrand 
in p.20p is norm continuous for z ^ M. If we set 

U "^^^ {z ^ X + iy\{x) < \y\<2{x)}, V '^^^ {z = x + iy\0 < \y\ < 2{x)}, 
then the norm of the integrand is dominated by 

c ^ |a™vl/(x)| — (a;)'"-2||ar|U^([i,2])li/(x + ty)+ 

m=0 

+c|5^+^*(x)|^|yr(^) ||r|U»([o,2])ly(:^ + *y)- 
Integrating with respect to y for > a + 1 yields the bound 

AT 

m — 1 



\\^{-Ao)\\L.(n)^Lnn)< I (^l^™*(^)K^) 



|a^+i'f(a;)|(a:)^) 



dx = lll'l 



JV+l- 



One may then prove that the operator ^i(~Ad), acting on L^'(il), is independent 
of TV > 1 + n/2 and of the cut-off function r in the definition of ^I^, see [12j. 
We now recall two lemma which will be useful when composing operators. 

Lemma 3.14 (Lemma 2.2.5, [T^). If 'i> E C^{M.) has support disjoint from the 
spectrum of —K^Ad then '^{—K^Ad) = 0. 
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Lemma 3.15 (Lemma 2.2.6, [12]). If^^i, *2 e C^(K); then {'^i^/2)i-h'^^D) = 

*l(-/i2Az5)*2(-/l2Az5). 



For the remaining part of the Appendix we prove the resolvent estimate p.22p 
from Lemma 13.111 If Rez > 0, this is nothing but a standard eUiptic estimate. The 
trouble comes with Rez < and getting close to the spectrum. In R", one may 
evaluate directly the convolution operator by proving its kernel to be in L^: this 
follows from 



N + ICn^ = sin^^V^(kl + ierf +cos2^V^(iei^-Nf, with z = \z\e^\ 



and a direct computation of LF' norms of d"{z + |^P)^^. By reflection, one then 
extends this estimate to the half-space case, with both Dirichlet and Neumann 
boundary conditions. By localizing estimates close to the boundary and flatten- 
ing, one may then obtain the desired estimate (|3.22p : such an approach is carried 
out in [T] in a greater generality (systems of Laplace equations, mixed boundary 
conditions), at the expense of fixing the angle 9 and not tracking explicit depen- 
dances on |2;| and 9. While (relatively) elementary, such a proof is, out of necessity, 
filled with lenghty calculations and most certainly does not provide the sharpest 
constant. It is worth noting, however, that it relies on standard elliptic techniques. 

To keep in line with the parabolic approach, we present a short proof, relying on 
the holomorphic nature of S{'w) in the half-plane Rew > 0. Remark that by our 
bound on S{t), t S IR.+ , the trivial bound on S{w), Kew > 0, and Stein's 
parameter version of complex interpolation, one may easily derive that S{w) is 
holomorphic in a sector around the positive real axis; but its angle will narrow 
with large or small p. However the argument may be refined and S^w) was proved 
to be holomorphic in the whole right half-plane in [12] , using in a crucial way the 
Gaussian nature of the heat kernel on domains ([10]). This was extented to more 
general settings in [5, where an explicit bound is stated: 



Then (|3.22p is a direct consequence of the following standard computation: recall 
the following formula, which is simply a Laplace transform. 



Jl 

where L can be chosen to be a half ray from the origin. Set z — re'^, w = pe*"^, 
then 



Now, if Rez > 0, we may take = and use estimates for the semi-group S{p). We 
would like to extend the range to the Rez < region, up to a thin sector around the 
negative real axis (|7r — 9\ < e); getting close to the spectrum is required if we want 
to define ^{-Ad) with ^ e Cg°{]0,+oo[). One picks (j) such that 216* -f 0| < tt, 
which ensures a decaying exponential in l3.25i provided we bound S{w) in L^. But 
the condition on cf) yields \4>\ < 7r/2, and the bound amounts to the holomorphy of 
S{w). The constant in p.24p translates into a (|2;|/|Imz|)" factor, while integration 
over p provides the remaining l/|Im2;| in p.22p . This concludes the proof. 



(3.24) 




(3.25) 
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